Drift, creep and pinning of a particle in a correlated 
random potential 



Heinz Horner 

Institut fiir Thcorctischc Physik, Universitat Heidelberg, Philosophenweg 19, D-69120 Heidelberg 
To be published in Z.Physik B. Submitted Aug. 1995. Revised Oct. 5. 1995. 



in 



o 
O 

> 
o 

00 

o 

ON 
1 



X3 
O 

o 



X 



Abstract. The motion of a particle in a correlated ran- 
dom potential under the influence of a driving force is 
investigated in mean field theory. The correlations of the 
disorder are characterized by a short distance cutoff and 
a power law decay with exponent 7 at large distances. 
Depending on temperature and 7 , drift with finite mo- 
bility, creep or pinning is found. This is in qualitative 
agreement with results in one dimension. This model is 
of interest not only in view of the motion of particles 
or manifolds in random media, it also improves the un- 
derstanding of glassy non-equilibrium dynamics in mean 
field models. The results, obtained by numerical inte- 
gration and analytic investigations of the various scaling 
regimes in this problem, are compared with previous pro- 
posals regarding the long time properties of such systems 
and with replica calculations. 



1 Introduction 

This paper deals with the motion of a particle in a corre- 
lated random potential under the influence of a driving 
force. The correlations of the disorder are character- 
ized by a short distance cutoff and a power law decay 
at larger distances with exponent 7. This problem has 
been solved in one dimension |l|, |[ |^ and, depending on 
7, drift with finite mobility, creep or pinning was found. 
Creep means that the particle moves with a mean ve- 
locity less than proportional to the driving force and in 
case of pinning the mean velocity is zero, unless the force 
exceeds some critical value. 

In N dimensions a mean field treatment becomes ex- 
act for — > 00 and there is a close formal similarity 
to spin glass problems with long range interactions. The 
situation without external driving force has been treated 
within replica theory j|J|] and also using a stochastic 
dynamics approach ^ Actually the more general 
case of the motion of a /^-dimensional manifold in an 
A — > cx) dimensional space was investigated. For finite 
A the mean field treatment is only an approximation. 
Nevertheless it might give some clue about a great va- 
riety of systems with disorder, for instance fiux lines in 
superconductors or interfaces in random field systems. 



The close formal relation to spin glass problems raises 
another set of questions. There has been a continuous 
interest in the dynamics of spin glasses and related sys- 
tems. The pioneering work of Sompolinsky and Zip- 
pelius ^, ^ concentrated in part on an understanding of 
the replica theory and the hierarchical replica symme- 
try breaking scheme proposed by Parisi ]lO| . The basic 
assumption was about the existence of a hierarchy of 
diverging time scales with ultrametric properties. This 
point has been investigated further |]ll| introducing a 
long time scale i associated with slow changes of the ran- 
dom interactions in this model. It could be shown that 
the above assumption meant that the long time contri- 
butions to the correlation and response functions have to 
be smooth functions oi x{t) = 1 — In i/ In f with i — > 00. 
It was later shown ||l2| that this led to inconsistencies 
which could not be resolved at that time. 

There exists a variety of models with disorder where 
a replica treatment requires only single step replica sym- 
metry breaking instead of the hierarchical scheme men- 
tioned above. The resulting phase transition can be con- 
sidered as discontinuous in contrast to the continuous 
transition observed in the SK-model and other cases re- 
quiring hierarchical replica symmetry breaking. A treat- 
ment via dynamics ||, ^ |l^, |l^ revealed a dynamic 
freezing temperature above the transition temperature 
obtained in replica theory. This is in contrast to the con- 
tinuous transitions where the same transition tempera- 
ture is found. Furthermore it could be shown |jl^ that 
the states which contribute to dynamics have a higher 
energy than the states relevant for the replica calcula- 
tion. Numerical simulations of a learning process in a 
perceptron with binary bonds [|4| , another system with 
discontinuous transition, indicated that the replica re- 
sult applies if the limit t ^ 00 is performed first and 
then the limit A ^ 00, whereas the results obtained via 
dynamics hold for the opposite order. This observation 
suggests that in the thermodynamic limit the particular 
states most relevant for the replica calculation cannot 
be reached within any finite time. This picture is con- 
sistent with the fact that finding the best groundstate 
in a spin glass or perfect learning in a perceptron with 
binary bonds is a combinatorial hard problem which can- 
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not be solved in polynomial time. The replica treatment, 
on the other hand, is purely static and does not refer to 
any kind of dynamics. For a given physical situation the 
appropriate order of limits has of course to be found, 
keeping in mind that the longest equilibration times are 
likely to diverge exponentially with some power of N 

This raises the question whether similar discrepan- 
cies also show up in systems with continuous transition. 
The problem of the motion of a particle in a random 
potential under the influence of a driving force is in sev- 
eral aspects an appropriate model system to address this 
question. Prescribing the drift velocity v and calculat- 
ing the force necessary to drive the particle, a long ex- 
ternal time scale is defined hy t ^ v~^. For times ex- 
ceeding i the system is supposed to reach a stationary 
non-equilibrium state, a so-called flow state, where cor- 
relation and response functions depend on differences of 
time only. This is of great advantage regarding numeri- 
cal as well as analytic work. Depending on the exponent 
7 and on temperature, a continuous as well as a discon- 
tinuous transition can be found and therefore both cases 
can be studied on a single model. 

The results presented in the following are based on 
numerical integrations of the dynamic mean field equa- 
tions of this model, extending over more than thirty or- 
ders of magnitude in time and velocity. Such a wide 
range is actually necessary in order to understand the 
different scaling and crossover regimes emerging in this 
problem. Such a wide range of time scales is not unreal- 
istic keeping in mind that phenomena might be observed 
on a scale of days or more with an intrinsic time scale 
of 10~^^ to 10~^^ seconds typical for vibrational or elec- 
tronic degrees of freedom. 

The numerical results are supplemented by exhaus- 
tive analytic work characterizing the different scaling 
regimes and focusing on their asymptotic properties and 
associated exponents. The characteristic time scales are 
obtained by matching between adjacent scaling regimes 
and are again governed by characteristic exponents. 

The key result of the present investigation is the ob- 
servation that the dynamics with a long but finite exter- 
nal time scale t is ruled by three regimes: 

i) The FDT-regime describing local equilibrium at short 

times. 

ii) The intermediate plateau regime for times t ^ 

with < ^ < 1. 

iii) The asymptotic regime for t ^ t. 

These regimes have additional substructures and associ- 
ated additional internal time scales. 

The plateau regime is the crucial link between short 
times and long times of the order of the external time 
scale. Its properties determine for instance whether 



creep or pinning is found or whether the longest inter- 
nal time scale is proportional to t or t ^+'' with 77 > 0. 
This question is of relevance for instance in the context 
of aging phenomena |l9j . 
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Figure 1: Phase diagram. O: Drift phase with v ^ F 
(driving force F, resulting velocity v). A: Creep phase 
with V - i^i/''. B and C: Pinning phases with v — Q 
for F < Fp. The phase diagram resulting from replica 
theory is also shown (doted lines) 

Depending on 7 and temperature, different phases 
are found. The phase diagram resulting from the present 
investigations is shown in Figj^. There is a phase with 
finite mobility (O), another one with creep behavior {A) 
and two more where pinning is observed {B and C). This 
is in qualitative agreement with the findings for the one- 
dimensional model which indicates that the transi- 
tions are not an artifact of the mean field treatment. In 
all phases, except the high temperature phase C, addi- 
tional intermediate time scales exist, which diverge for 
V 0. There is, however, no ultrametric organization of 
time scales. This is in contrast to the common proposal 
0, ^ |l5). Moreover, the phase diagram suggested by 
replica theory |4[ ^ or dynamics within the hypothesis 
of ultrametric time scales ||, Q differs from the one ob- 
tained in the present investigation. For instance, 7=1 
has been predicted for the A-B phase boundary in con- 
trast to 7 = 1.3044 found at present. This leads to the 
conclusion that the states dominating the replica theory 
with broken symmetry are again not the same as those 
relevant for dynamics on long but finite time scales. 

As a consequence of the assumption of ultrametric- 
ity in the organization of diverging time scales, the 
resulting correlation and response functions were not 
unique |ll|. The present investigation yields on 
the other hand unique results and this problem is 
therefore removed. The inconsistencies [|2| mentioned 
above are also resolved in the sense that the proposed 
reparametrization of time x{t) = 1 — Int/lnt is not ap- 
propriate. 

There has recently been a growing interest in 
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the dynamics of glassy systems evolving from a non- 
equilibrium initial state |l|, |2|, |2|, |24|, |2|, |2|l . The 
situation there is more complicated because the problem 
is no longer homogeneous in time and response and cor- 
relation functions depend on two time arguments t and 
t' . If, however, t — t' ^ t', one expects a behavior similar 
to the one investigated at present with t' taking over the 
role of i. Especially if both t — t' and t' diverge with 
{t — t')/t' 0, one should also expect the appearance of 
intermediate time scales and associated scaling regimes. 
This possibility has not been considered at full depth in 
the above work. 

For short time t — t' correlation and response func- 
tions are related by fluctuation-dissipation theorems 
(FDT) . This is also the case in the present investigation. 
The different phases have, however, different scenarios of 
how FDT's are violated at longer time scales. Since this 
happens for times 1 <^ t <^ t, similar phases are ex- 
pected for the above mentioned relaxation processes. In 
this context the non-equilibrium dynamics of the spher- 
ical SK-model investigated recently by Cugliandolo and 
Dean is of great interest, because this model does 
not require replica symmetry breaking in its low tem- 
perature phase. 

The present paper is organized as follows: Section 
2 contains the definition of the model and the dynamic 
mean field equations, which have the form of coupled 
nonlinear integro-differential equations. It also contains 
the definition of effective time dependent exponents, 
which are widely used in the following. Section 3 is 
devoted to the numerical integration of the dynamic 
mean field equations, to the presentation of results and 
to a preliminary identification of different time scales 
and scaling regimes. Section 4 contains analytic results 
and estimates regarding time scales and scaling prop- 
erties. It starts with the FDT-solution valid for short 
time and continues with a discussion of various convolu- 
tion type integrals entering the self-energies. In bypass- 
ing the QFDT-solution |^ and the hierarchical solution 
are brought up. It is then shown how to evaluate 
the above integrals using the time dependent effective 
exponents. This leads to a reformulation of the original 
mean field equations in terms of a set of 15 coupled ordi- 
nary and first order differential equations. This appears 
complicated but it allows a discussion of different scaling 
regimes in which only closed subsets of these equations 
are relevant. These regimes are in the order of increas- 
ing time the FDT-regime, the plateau regime, where the 
correlation function q(t) stays close to the asymptotic 
value qc of the FDT-solution, and finally the asymptotic 
regime for times of the order of the external time scale 
v~^. A concluding discussion follows in Section 5 and a 
brief derivation of the dynamic mean field equations is 
given in an appendix. 



2 Formulation of the problem 
and dynamic mean field theory 

2.1 The model 

A particle is investigated which moves under the influ- 
ence of an external force F = {^/NF, 0, • • • , 0} in a ran- 
dom potential V{q). The coordinates of the particle 
are written as g = {gi + ^/Nvt, P2, ■ • • , gjv} assuming a 
mean velocity v = {v^w, 0, • • • , O}. This means that 
the components Qi are measured in a frame moving with 
velocity ^/Nv along the 1-direction. The scaling with 
\/N ensures a nontrivial limit N ^ oo. The system is 
at a temperature T = {3^^ and its Hamiltonian is 

H^V{Q)+^^^i^Q'-F■Q, (2.1) 

where a confining potential ^^jlqQ^ might be added to the 
random potential and the potential of the driving force. 
The following investigation is, however, restricted pri- 
marily to /io = 0. The motion of the particle is governed 
by a Langevin equation 

dtQc. = -P^+^<. (2.2) 

OQa 

with Gaussian white noise 

{Ui)ip{t')) = 'i5^p5{t-t'). (2.3) 

The quenched disorder is also assumed to be Gaussian 
with 



v{q) = 

v{Q)V{e') = ~Nf{j^{e-e'f) (2.4) 

and 

f{x)^,^^^^{l + xf-'^ . (2.5) 

The strength and the short range cutoff of the disorder as 
well as the diffusion constant are set to one. This model 
differs from the one studied earlier ^, ^, Q only by 
the addition of the driving force. 

2.2 Dynamic mean field theory 

The order parameters of the dynamic mean field theory 
are the correlation function 



OL 

and the response function 



r{t-t')^-Y^5{g^[t))/5F^{t') . (2.7) 
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They obey the mean field equations |^ 

dtr{t) = —fir{t) + / ds ■w{s)r{s)r{t — s) 
Jo 

and 

dtq{t) ^2- fiq{t)+ I dsw{s)r{s)q{t - s) 
Jo 

- ds^2\W{t + s)~W{s)\r[s) 

— [w(t + s)r(f + s) — w(s)r(s)] (7(s)| 

W{t) = 2l3^f{q{t)+vh^) 



(2.8) 



with 



and 



M = /^o + 



ds w{s)r{s) . 



(2.9) 

(2.10) 
(2.11) 

(2.12) 



A brief derivation is given in the appendix. In the fol- 
lowing /xo = is assumed. 

In thermal equilibrium correlation and response func- 
tions are related by a fluctuation-dissipation theorem 
(FDT) which reads dtq{t) = 2r{t). This is violated in a 
non-equilibrium situation and 



nit) 



1 



2r{t) 



dtq{t) - 1 



(2.13) 



is introduced as measure of this FDT-violation. 



Using this in ( |2.10D and ( |2.1l| ) 

dtW{t) = -{[1 + n{t)]rit) + v^t^wit) (2.14) 



is found. Combining (|2.8|), (2.9) and (2.13) one obtains 

rt 



r{t)dtn{t) 



dsw{s)r{s)r{t — s) 
X [n{t) — n{t — s)] 







ds- 



^w{t + s)r(t + s)r{s) 
X [n{t + s) — n(s)] 

+ v^{t + s)w{t + s)r{s)^ (2.15) 



and this equation can now be used instead of one of the 
original mean field equations (2.8) or (2.9). 

Prescribing the drift velocity v the average of the 
force necessary to drive the particle is 



ds 



(2.16) 



2.3 Effective exponents 

It is convenient to define 



and then ( 2.13 ) is written as 

tdtq{t) = 2{l + n{t)}g{t). 



(2.17) 



(2.18) 



In order to present the results of the numerical in- 
tegration of the mean field equations effective time de- 
pendent exponents are introduced. They are also used 
in the evaluation of the convolution type integrals in th e 
self-energies of the mean field equations (2.S) and ( 2.15 ). 
The first two exponents are 



and 
With 



z/(t) = tdt\ng{t) 
a{t) = -tdt\nw(t) . 



k{t) = tdtn(t) 



the third exponent 



K{t) = tdt\nk{t) 



(2.19) 

(2.20) 
(2.21) 

(2.22) 



is defined. 



3 Numerical results 

3.1 Numerical procedure 

The dynamic mean field equations to be solved consist of 
two integro-differential equations (|2.8|)and (2.15), a first 



which is derived in the appendix. 



order differential equation ( ^.13 ), the definition oiw{t) 
( 2.10| ) and (2.5), and the integral (2.12) determining /x. 
The second integral in (2.15) and the integral ( 2.12] ) ex- 
tend overall times and therefore, evaluating the solution 
at some t requires the knowledge of the solution at all 
times, not only a± s < t. This means that the equa- 
tions have to be iterated. This is the price one has to 
pay for investigating a steady state situation. A study 
of the relaxation from a non-equilibrium initial state on 
the other hand is free of this problem, but it requires to 
deal with functions depending on two time arguments. 

The wide span of time arguments ranging from 
10""* to 10^^ requires a nonuniform discretization. The 
present calculation uses a homogeneous discretization of 
hit. The convolution type integrals in (2.8) and (2.15) 
can be evaluated by numerical integration using the val- 
ues stored at the points of the above grid for one of the 
functions of the integrand and an interpolation for the 
other function. Alternatively, an approximative evalua- 
tion of the integrals using effective exponents, described 
in Section 4.4 is possible. For the actual calculation both 
methods have been combined. 
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Figure 2: Correlation function q{t) in phase A 



Figure 4: •n{t) in phase A 
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Figure 3: Correlation function q{t) in phase B 



Figure 5: n{t) in phase B 



The mean field equations in the form written in Sec- 
tion 4 consist of several coupled nonlinear differential 
equations of first order. They can be integrated forward 
or backward in time. In addition they have to be iter- 
ated. The situation resembles in some sense the problem 
of solving 

Ax [t) 



dt 



(3.1) 



Integrating forward in t the solution typically ap- 
proaches a fixed point f{x) = with f'{x) < 0, whereas 
a fixed point with /'(x) > is reached integrating back- 
ward in t. Selecting an appropriate initial value a;(to) 
somewhere between a fixed point with f'{x) > and an 
adjacent one with f'{x) < 0, the solution can easily be 
found by integration in both directions. For the com- 
plete set of equations such fixed point situations show 
up at the crossover between the various scaling regimes 
discussed below. 

3.2 Results 

The following figures show selected results obtained for 
{7 = 1.5; T = 0.4} and {7 = 0.75; T = 0.7}, respec- 
tively. The first set of variables corresponds to a situ- 



ation where the replica calculation ||] requires single 
step replica symmetry breaking. This corresponds to a 
point in phase A of the phase diagram, Figj^. The sec- 
ond set corresponds to a point in phase B and the cor- 
responding replica treatment requires continuous replica 
symmetry breaking. The following figures show results 
obtained for various values of the drift velocity ranging 
from V ~ 10~^ to V — 10^'^". In the following various 
velocity dependent characteristic time scales are intro- 
duced. Their values for v — 10~'^° are marked in the 
figures. 

Figsjl and || show the cor rela tion function q{t). The 
plateau value Qc is given by ( |4.3| ) . In view of the argu- 
ment q{t) + v'^t^ in ( I2.IOD the value of vH"^ for v = 10-^° 
is also shown. 

The function n{t) defined in ( 2.13 ) is shown in Figs.^ 
and p. It indicates the violation of the FDT, which holds 
for t < tp. It develops a first plateau value for tp < t < ta 
and a second one for t > t'^. In phase A this second value 
obviously depends on v. In phase B the first plateau is 
not yet fully developed, even for v = 10""^°. 

The time dependent effective exponents i^{t), ( 2.19|) , 
k(<), ( |2.22| ), and a{t), ( |2.20| ), are shown in Figs .| to |ll|. 
Their properties are discussed below. 
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Figure 6: Effective exponent v{t) in phase A 
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Figure 8: Effective exponent K{t) in phase A 
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Figure 7: Effective exponent i/{t) in phase B 

From the data several characteristic time scales can 

as 



be extracted. They are shown in Figs.|l^ and 13 
functions of the velocity v. The time tp is defined by 
q{tp) = Qc and it characterizes the center of the plateau. 
The next time scale is relevant only for phase B and 
it is given by Kit^) = 0. The time scale ta is defined 
by 9(^0) = 2gc and ijj by q{t'^) = 2v'^t'^2. A power law 
dependence is found for v ^ 0. 

3.3 Scaling regimes 

The numerical results and analytic considerations de- 
scribed in Section 4 reveal the existence of various scal- 
ing regimes in the limit z) — > 0. 

3.3.1 FDT-regime 

For finite t ^ to one finds n{t) ^ 1 and fluctuation 
dissipation theorems hold. This can be understood as 
a situation where the particle stays within one valley of 
the energy landscape. The correlation function and n{t) 
can be written as 

q{t) = q{t/to) 
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Figure 9: Effective exponent K{t) in phase B 

n{t) ^b{v)n{t/to) (3.2) 

with b{v) — > for V — > 0. The time scale to and the 
functions q{x) and n{x) do not depend on v. This regime 
is referred to as the FDT-regime. 

3.3.2 Plateau-regime 

With increasing t a plateau regime is found where q{t) ~ 
qc- The corresponding time scale is tp = tp{v) defined 
by q{tp) = qc and the above functions obey the scaling 
form 



qit) = gc + a{v) q{t/tp) 
n{t) = ric + n{tltp) 



(3.3) 



which is consistent with v{t) = i){t/tp) and K{t) = 
k{t/tp). 

The function n(i), which is a measure of FDT- viola- 
tion, reaches a plateau value ric for t ^ tp. In phase A 
this plateau extends beyond the limits of the g-plateau. 

In phase B the plateau of n{t) is not yet very pro- 
nounced, even at u = 10~^°, but it can be seen that the 
upper limits of the n-plateau and the g-plateau essen- 
tially coincide. The center of the n-plateau defines an 
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Figure 10: Effective exponent a{t) in phase A Figure 12: Characteristic time scales in phase A 




Figure 11: Effective exponent a{t) in phase B 



Figure 13: Characteristic time scales in phase B 



additional time scale tx = tx{v) where n(tx) — tic or 
K{tx) — 0. For t ^ tx a. scaling form 

n{t) ^ nc + c{v)nx{c{v) Int/tx)) (3.4) 

will be derived. 

In phase C n{t) <^ 1 holds up to the upper end of 
the (/-plateau and consequently the FDT-solution holds 
in the whole plateau regime. 

3.3.3 Asymptotic regime 

The lower end of the asymptotic regime is marked by a 
time scale ta — ta{v) where q(ta) ~ Qc ^ Ic- An appro- 
priate choice is ta such that q{ta) = 2qc. For t ^ ta & 
scaling form 



q{t) = q{t/ta) 

n(t) ^ nc + b{v)n{t/ta) 



(3.5) 



with b{v) = 1 in phase B and C is found. 

Another time scale t'a — t'a{v) can be defined such 
that q{t'a) — (wt'a) . In phase B and C both are propor- 



tional to 



whereas ta{v) ^ v ^^"^ and t'a{v) ^ v ^ ^' 



with rj > is found in phase A. This requires for t ^ t'a 
the scaling form 



n{t)=nc + b'{v)n'it/t'a). 



(3.6) 



The existence of the exponent rj > is connected to the 
creep behavior observed in phase A. 

The above results now have to be verified by analytic 
investigations. This is done in the following sections. 

4 Analytic results 



is tp with 



4.1 FDT-solution 

The shortest w-dependent time scale for v 
1 ^ <C For t <^ tp one expects a solution which 
obeys the FDT and therefore n{t) = 0. This solution 
holds in phase O for v = and describes equilibrium 
within a single valley of the energy landscape in the other 
phases. The discussio n o f t his so lution follo ws standard 
arguments @. With (U), ( |2.14| ) and (^.18| ) one obtains 



dtq{t) = 2 - ^lq{t) 



+ / ds w{s)r{s)q{t — s) , 
Jo 



(4.1) 
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because the second integral in (2.9) vanishes. In order 
to demonstrate this, one realizes that its main contribu- 
tions come from s ^ t where n{t) sa 0. With ( p. 17 ) 
and dOI ) dtq{t) « 2r{t) and with ( |2l^ ) dtW{t) « 
—w{t)r{t). This allows to write the integrand as a com- 
plete derivative with respect to s and to evaluate the 
integral resulting in 0, since q(0) = and W{t) for 
t — > oo. 



Eqs.(2_^) and (4J) cannot be solved in closed form. 
It is, however, sufficient to investigate the solution for 
t tp, where q{t) qc, dtq{t) ~^ and dtr{t) 0. 

In leading order the integrals in both equations can 
be evaluated by taking into account only the contribu- 
tions near the upper and lower bound, respectively. This 
yields from ( ^ ) 



2/3M/'(0)-/'(gc)-9c/"(<?c)} (4.2) 



and from (4.1) 



nq,)ql = -T^ 



(4.3) 



The second e qua tion is used t o determine qc ■ E li mina- 
tion of fi in using (|2.12D yields with ( ^.10[ ), (l2.1l|) 

and (|2T^) 



qcf'iqc) = - 



2/32 



ds w{s) r(s) 



This condition involving the solutions at long time scales 
will be used later. 

In order to give estimates of the corrections to the 
leading order one assumes 



tr{t)^g{t)~^-ii^oqc (t/to)" 



(4.5) 



with uq < 0. The time scale to has yet to be determined. 
This leads to 



lit) 



ds r(s) 



,{l^{t/tor}. (4.6) 



This is now used in an analysis of the corrections in 
(2.8) where the leading order contains terms ^ f^"'^. 



The next to leading order ~ t'^'^" ^ results in 



f"{qc) I ds {ris)rit ~ s) - 2rit)ris)} 

{f'iqc) + if"'{qc)qcr{t)}{qe - qit)} . (4.7) 



With ( [4. 51) the integral on the left hand side can be 
evaluated. Introducing 



qcfiqc) 
2/"(9c) 



(4.8) 



one obtains 



r2(i + 1^0) 



— Rn 



(4.9) 



r(i + 2i.o) 

which allows to determine the exponent vq. 

In order to get an estimate of the time scale to one 
can use an interpolation formula 

q{t)^qc[l-{l + t/tor} (4.10) 

which gives the correct asymptotic behavior for t tp. 
The requirement q{t) 2t for t — > yields 

to~-^iyoqc- (4.11) 

In the ergodic high temperature phase O fo r v = 
the FDT holds for all times and therefore with ( 2.12 ) 

(4.12) 



M = 



2^2 



At the transition to one of the nonergodic phases the 
plateau develops and one finds with (4.2) and (4.3) 



T! - qj'iqc) , 
f'iqc) = -qcfiqc) • 



(4.13) 



A solution of these equations with finite Tc exists only 
for 7 > 1. It is 

1 



qc = 



7-1 



-7/2 



(7-1) 



(7-l)/2 



(4.14) 



This is shown in Fig.|l|. The corresponding value ob- 
tained from replica theory |j 



clRSB 



67 



(4.15) 



is always below the above value. 

For T > Tc in phase O the FDT holds for all times 
and the force (2.16) is 



fiF = v^l + dtWit)^ 



(4.16) 



For t ^ 00 (2.8) yie lds r {t) — > Too and qit) ^ t. There- 
fore the integral in (4.16) converges for 7 > 1 and there 
exists a finite friction constant for v 

4.2 Integrals and counterterms 

For t 3> 1 the leading contributions to the integrals in 
( |2.8| ) and (2. If;) come from regions near the boundaries 
of integration. Subtracting those leads to residual in- 
tegrals and a reformulation of the dynamic mean field 
equations. The following integrals are introduced: 



J{t) = ^ ds |u;(s)r(s)r(i - s) - w(t)r(t)r(t - s) 

-wis)ris)rit) + swis)ris)dfrit)^ , (4.17) 
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K{t)^ / dsw(t + s)r{t + s)r{s) 
Jo 

X |n(s) — n(s + t)| 

;u'(t — s)r{t — s)^r{s)n{s) 
-{s)n{t) + (t~ s)r{t)dtn{t)^ (4.18) 



ds 'i 



and 

/•OO 

[/(i) / ds + s)w{t + s)r(s) 

Jo 

Furthermore it is appropriate to define 



(4.19) 



/OO 
ds w{s)r{s) - w{t) J dsr{s) (4.20) 



and 



Z{t)^l+ dssw{s)r{s). (4.21) 

"'0 



This aUows to rewrite the mean field equation (2 
in the form 



Z{t) dtr{t) ^ J{t) - D{t)r{t) (4.22) 



and ( |2.15| ) as 



Z{t)r{t)dtn(t) = Ult)-K{t). 



(4.23) 



4.3 QFDT- and hierarchical solution 



With ( [4.20D and n(i) w for t < the condition (^ 
means D{tp) — 0. 

Depending on 7 for t; = a QFDT-solution and a 
hierarchical solution, respectively, have been proposed 
1^, |7j. These solutions can easily be obtained from the 
above equations. Especially the second scheme requires, 
however, assumptions which are not fulfilled. Neverthe- 
less they are presented here for further reference. 

4.3.1 QFDT-solution 

In the QF DT-solution n{t) — uq for t > tp is assumed. 
From (44) one finds 



nQ + l = - 



file) 

qcfiqc] 



(4.24) 



Without further analysis of the regime t > tp nothing 
can be said about the range of validity of this solution 
or about the force F. This solution was proposed to 
be valid for 7 > 1 and it shows some similarity to the 
IRSB-calculation g |. 



4.3.2 Hierarchical solution 



Assume that the integral in (2_^) for i > is com 



pletely determined by its contributions from the up- 
per and lower bound, respectively. This means that 
J{t) « 0. As it turns out this is the essence of the pro- 
posal of the existence of an ultrametric organization of 
long time scales. Neglecting the derivative with respect 
to t in ( |4.22D results in D{t) fa for t > tp. 
Defining 



(4.25) 



Eq.(g^ yields 



dq' nq')m{q') = f"{q) / dq' m{q') . (4.26) 



(4.27) 



This leads to the differential equation 

/""(<?) 3/"'(g) 



which is solved for q > qc 

1 + . 



771(9) = ™('7c) 



1 



2/"(g) 



-i(i-7) 



(4.28) 



using (|2J). 

The integration constant m(qc 
(|4.26|) with q^qc 



is obtained from 



m{qc 



f"'iqc)qc 
2/"(gc) 



R„ 



(4.29) 



with R„i given in (4.8). This is the solution obtained 
previously, assuming a hierarchical structure of long time 
scales 0. 

4.4 Evaluation of integrals via effective 
exponents 

The main problem in a discussion of the solutions of 
the mean field equations and also in their numerical in- 
tegration is in the evaluation of the integrals listed in 
the previous section. In the following an approximative 
scheme is proposed, which turns out to be very accurate 
over the whole range of t. This scheme is based on the 
effective time dependent exponents, which have been in- 
troduced in Section 2.3. It allows to write and evaluate 
the integrals in the form 



dx ; 



r(a)r(6) 

T{a + b) 



(4.30) 



As a first step in this program one introduces " dimen- 
sionless" quantities instead of the original ones given in 
(p7|) to (E2l|): 



J{t) = 



tJ{t) 
w{t)9Ht) 



(4.31) 
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and 



m 

m 



tK{t) 



tUjt) 
k{t)w{t)g^t) ' 

m 



w{t)g{t) 



(4.32) 
(4.33) 
(4.34) 

(4.35) 



twit)g{t) ■ 

The functions in the integrands are approximated by 



sr{s)=g{s)^g{t)(^^y 



nj{s) « wit) (0 



-ait) 



and 



n{s) w n{t) 



m 



1- 



(4.36) 



(4.37) 



(4.38) 



This ansatz fulfills ( [ITgl ), (|^, ( ^^ and (|^ 
for s = t. The integrals J(t) and K{t) are no w of t he 
form proposed above (for the first integral in ( [f.lS| ) a 
trivial substituti on of variables is necessary) and can be 
evaluated using dOOl ) provided a > and 6 > 0. Oth- 
erwise counterterms are required cancelling the poles of 



the F-functions. The last three terms in (4.17) are of this 
kind. Depending on the values of the effective exponents 
further counterterms might be required. 

This yields (in the following the arguments t are 
dropped) 



J'. 



1 1 i^-l 



v V — a. v + 1 ~ a 



(4.39) 



and 



X ^ 1 J r(i/ + K)r(l + a-2iy~K) 



K I r(l + a — ly) 

r(i/)r(i + a - - k) 

r(l + a — h' — k) 

r(i/-a)r(i/) r(j/ - a)r(t/ + k) 

r{2v-a) r{2iy + K-a) 



V — a + 1 



(4.40) 



by 



The integral (4.33) is given for > and a — v > 1 

-^^r(.)r(a-.-i) ^^^^^^ 



otherwise it is approximated by 

v'^tZ{oo) 



U : 



kw g"^ 



For a — u > I ( 4.35 ) becomes 

~ twg 

otherwise 



1 



(4.42) 



(4.43) 



(4.44) 



The mean field equations (4.22) and (4.23) are rewrit- 
ten as 

{iy-l}Z=J-D (4.45) 

and 

Z = U-k. (4.46) 



The exponent a, ( 2.20 ), is 



-2g{\ + n + u] 



f"'iq + vH^) 
f"{q + vH"^) 



with 



2 j-2 



Differentiation of ( [1.20D yields with ( |4.34| ) 
tdtb = {a - v]b - 2{l - [1 + n + u] R] 
with 

R ^ - •'„„" ' I dsr{s). 



(4.47) 



(4. 



(4.49) 



(4.50) 



kg 



T{a-1) 



The complete set of mean field equations now in- 
volves 15 functions of time: q, r, g, n, k, v, k, a, u, 
R, D, J, K , U an d Z. T he corr es pondi ng 15 e qu ation s 
are: (M, ( ^ ^ (B , (HI), S), (H), 

(|l4l|/ p^ ) (fl.43| / [4.44| ), ( [4.45D , ( |4.46D , ( Pf ), glj ), 
(|l49| ) and (|4.50|). 

Some of the variables can easily be eliminated, but 
even then one is left with 5 differential equations of first 
order and 2 implicit ordinary equations for 7 of the above 
functions of time. This certainly looks complicated. On 
the other hand, in each of the scaling regimes to be dis- 
cussed below only a subset of variables and equations 
has to be looked at. 



4.5 The plateau regime 

Having discussed the FDT-solution in Section 4.1 al- 
ready, the investigation of the plateau regime follows 
next. This regime is characterized hy \q{t) — qc\ <^ qc- 
Further simplifications are due to t^v ^ qc- Conse- 
quently R defined in ( 4.50D can be replaced by the con- 
stant Rm given in (4.8). Furthermore {/ w and u « 
can be used. From (2.20) with (2.10) one finds a w 
and then ( |4.44| ) yields Z w jj^. 
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-1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 14: i^(-D) and —Di/{D) as functions of D, see 
text. 



Eqs.( [4.39D , ( |4.45| ) and ( [4. 44]) determine v = v{D) as 
a functi on of Z), sh own i n Fig.|l4|. 

Eqs.( [4.40|) and ( [4.46[ ) are solved by 



= 1 - ■iv{t) . 

Eq. ( 4.49D now reads 

tdtD{t) = ^D{b{t))D{t) 

-2{1- [l + n(<)]i?„,} 



(4.51) 



(4.52) 



The remaining set of equations to be solved consists 
of ( p^ , ( |l2^ ), ( |44^ ), ( p5| ), ( ^Hl ) and ( |I52| ). 

Time enters only via derivatives of the form tdt and 
therefore a family of scale invariant solutions exists. This 
means that the whole family can be written as a set of 
functions depending only on t/tp with arbitrary tp. This 
is in accordance with the scaling form proposed in (3.3). 

The time scale tp 3> can now be fixed such that 
D{tp) = and as a c onsequence v{tp) — and K{tp) — 
1. Integrating ( 4.52 ) down to t ^ a stable fixpoint 
D{t) Do > is reached with 



HDojDo ^ 2{1 - R^} 



(4.53) 



and n{t) ^ 0. Eq.(gJ) yields iy{Da) = uq. 

For to t ^ tp the function q{t) can be matched to 
the asymptotic FDT-solution (4.6) and 



d{v) = {tp/toY 



(4.54) 



results for the scale factor in ( p.3[ ) where tp — tp{v) has 
to be determined later by matching to the asymptotic 
regime. It can be shown that the above fixing of tp is 
equivalent to the original definition q(tp) — within a 
factor close to 1. 

The discussion of the asymptotes for t tp, still 
within \q{t) — (7c| <C qc, is more inv olved . Neglecting the 
derivative with respect to time in ( 4.52 ) one gets 



-Dv{D) = 2{l-(l+n)R„,} 



(4.55) 



which now determines D = Di{n), vi ~ v{Di[nYj and 
Ki = k{bi{n)) = 1 - iv{bi{n)). 

There is a special value rix such that ki^Diijix)) — 
with bi{n^) = b^ Ri -0.700084. This value is found 
from ([4.52D as 



1 



ID, 



R„ 



1, 



(4.56) 



Expanding 

R{bi{n)) K {n - n^) K (4.57) 
with k' > 0, Eqs.( |2.2i ) and ( |2.22 ) can be integrated in 



closed form. The real solutions are 



and 



2 c 

i[t) coth (c \n{t/t,)) 



n(t) = Tlx ~\ tan (c \n{t/tx)) 



(4.58) 



(4.59) 



with constants of integration c and tx- 

There are now t wo ch aracteristic values for n{ t), th e 
value uq defined in ( 4.24 ) and tlx defined above in ( 4.56| ). 
Both are functions of T and 7. 

The relative magnitude of ng and determines 
which of the above solutions has to be selected, as will be 
discussed in the following. This, on the other hand, de- 
termines ultimately whether creep or pinning is observed 
and therefore the different phases are characterized by 
their values of nq and Ux- Phase A is defined as the 
region where > nq > 0, phase B has nq > Ux > Q 
and nq > > n, holds in phase C. The ergodic phase O 
with finite friction has already been discussed in Section 
4.1 

4.5.1 Phase A 

The numerical results shown in Fig.^ indicate a plateau 
of n{t) extending to times beyond the upper boundary 
of the g-plateau. This suggests that the solution (4.58) 



has to be used for t > tp. The choice tx ~ tp yields the 



scaling form (3.3). For t ^ tp the requirement n{t) 
< i^x determines 



K 



Hi 

2 



(4.60) 



with Ki — k{bi{nq)^ . The parameter nq is later shown 
to ag ree with the value obtained for the QFDT-solution 
(HI)- 

For t ^ tp the above equations yield with vi = 
i>(bi{nq)) > and ki = 1 - 3i/i < 



q{t) qc + {tp/toY" {t/tpY^ qi 
n{t) nq - {t/tpY^ hi . 



(4.61) 



For w — )■ the constants qi and hi depend only on 7 and 
T. 



11 



For t < U 



n{t) 



no 



(t/tp) 



(4.62) 



with constant fiQ and kq > 0. Matching at to 'ti t <ti tp 
yields the scale factor 



b{v) = {to/tpY 



(4.63) 



in (3^). The time scale tp — tp(v) is still open and has 
to be determined later by matching to the asymptotic 
regime. The same holds for ng. 

4.5.2 Phase B 

The plateau of n{t) in this phase extends only to the 
upper boundary of the g-plateau, as shown in Fig.^. This 
indicates that solution (4.59) now has to be used for 
t > tp. Its range of validity is restricted to — ^tt < 
cln(t/t^) < iTT or e~''/'^H^ < t < e^/'^H^. Later it is 
shown that this range increases for u — s- and therefore 
c = c{v) — > for u — > 0. For t ^ tp the solution should 
not depend on c, which is the case for the choice 



,7r/2c 



Up . 



(4.64) 



For tp <^t <^ e^l'^ tp the second term in ( 4.59 ) is a small 
correction and the proposed scaling form (|3.3) and (3.4) 
holds over the whole range to <^t e^^'^ tp. 

With Ki — K(^Di{nx)) = and v \ — v {D\[nx]) — \, 
q{t) has again the asymptotic fo rm (I61I) for t ^ tp and 
the scale factor b{v) is given by ( 4.63| ). In this phase the 
quantities to be determined later by matching to the 
asymptotic regime are tp — tp{v) and t^ — tx{v). 

4.5.3 Phase C 



In this phase with n^; < again solution (4.58) has to 
be used. Since now ki — k(^Di{nx)) > and n{tp) ^ 1, 
the appropriate choice of the parameters is 



2c 



This yields for t^ ^ tp and tp <^ t <^ t^ 
n{t) « -2n, {t/t^y < 1 



(4.65) 



(4.66) 



which means that (^^) is fulfilled with n{t) « 0. 

The asymptotic form of q{t) is again (4.61) with vi = 
i'(-Di(0)). Again the quantities to be determined later 
by matching to the asymptotic regime are tp = tp{v) and 

4.6 The asymptotic regime 

The discussion so far did not depend on the actual choice 
of the drift velocity v. On the other hand, there are 
parameters not determined yet. These are the time scale 
tp{v) for all phases and the second time scale tx{v) for 



phase B and C. For phase A it has to be shown that ng 
actually agrees with the value obtained in the QFDT- 



solution, (4.24). 

The velocit y v enters t he full set of mean field equa- 
tions in U{t), ( 4.41 / 4.42 ), at various places in the ar- 
gument of f(q{t) + v'^t^ ) and its derivatives, and in the 
definition of u{t), ( 4.48 ). In order to obtain the scaling 
form proposed in ( |3.5| ) for t ^ ta it is necessary that (7 (t) 
can be written as U{t) — U{t/ta), which is the case for 



h{v) 



= 1 



(4.67) 



assuming the proposed scaling form in ( 4.41 / 4.42 ) . 

For t <^ ta the solution in the asymptotic regime 
has to match the solution valid in the plateau regime 
obtained in the preceding section for t^ tp. This means 



^ <lc 



(4.68) 



with constant go- Comparison with ( 4.61 ) yields 

(tp/toY" t;--^ - t-'^' ; tp-tl^ti (4-69) 

with 



c 



1^1 



(4.70) 



I'l - vo 

The discussion of n{t) has to be done for the different 
phases separately. 



4.6.1 Phase A 

Following ( |3.5| ) one can write 

n{t) HQ ~ {t/taY^ brio 



(4.71) 



with constant no and b — b{v). Comparison with (4.61) 
results in 

b ~ {ta/tpr^ . (4.72) 



Using this in ( 4.67 ) one finds with ( 4.69 ) 



with 



ta 

tp 

V 



-i+n ,n 

Ur 





^_(i_,)C^i-(i-,)C 



Vo Kl 



(4.73) 
(4.74) 



2(i^i - Vo) + t^oKi ' 

Since i/i > 0, i/q < and ki < one finds 77 > 0. 

For t ^ ta the remaining w-dependence in the mean 
field equations is not of relevance, since 



and 



^{vtof\tltaf <^q{t/ta) (4.75) 



u{t) = {v toY\t/taf/-g{t/ta) « 1 (4.76) 



with uit) defined in (4.48) and g(t/ta) = git) 
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For t ^ ta one obtains i'{t) 1, r{t) fi/ta and 
q(t) — > qi t/ta where fi and qi are constants. This fol- 
lows immediately from (4.43) and (4.45) realizing that 
Z(t) diverges ~ whereas the right hand side of 

( 4.45 ) remains finite for t ta- This defines an even 
longer time scale 

t'a-v-'-^t-^ (4.77) 



where q{t'^) = 

For ta < t < t' Eq. ( ^.20| ) yields a 7 + 1. In- 
vestigating Eq.( |4.46| ) one obser ves that Z{t), ( [4.43D , di- 



verges ~ t"-^ and U{t), ( |4.41| ), - t^-", whereas K{t), 
( 4.40| ) remains finite. This means K{t) « 3 — a(t) and 
K(t) K2 = 2 — 7. For 7 < 2 the exponent K(f) increases 
with increasing t reaching a value K2 > and therefore 
n{t) also starts to increase again for ta ^ t ^ t'^ reach- 
ing a new constant value 712 (w) for t 3> tj,. Matching 
with (|4.7lD and ([4.72[) yields 



n(t) 



,2),(1-K2) 



(4.78) 



and therefore n2{v) — 

2 this second plateau of n{t) is missing, the conclusions 
are, however, unchanged. 

In order to give an estimate of the driving force one 
has to investigate the behavior of a{t), Eq.( 2.20| ), first. 
In the plateau region a{tp) « was found. Around t ^ 
ta it starts to increase and reaches a value a{t) « 7 + 1 
for ta t <^ t'a- Around t ^ t'a it starts to increase 
again reaching its asymptotic value a{t) — 2(7 4-1) for 
t ^ t'a- The integral in the expression ( 2.16 ) for the 
driving force gets its main contribution from s such that 
1 -I- i/(s) — a{s) = 0. For 7 > 1 this is fulfilled for 
s ^ ta and n{t) is well approximated by ng. This yields 
the following relation between mean velocity and driving 
force: 

F^ivtoYF (4.79) 



with 



F = 



1 + J^Q Jo 



dx X f" {q[x)) d^q{x) (4.80) 



which does not depend on v. For given force the velocity 
increases slowly with increasing force according to 



(4.81) 



This behavior is a form of creep. 

In the one-dimensional case creep has also been 
found for i < 7 < 1 with 



(4.82) 



with fi = 2(7 — l)/(27 — 1). A power law dependence of 
the above form (4.81) is obtained for 7=1. Note that a 



different definition of 7 is used in |^ . The results quoted 
here refer to the present definition (2.5). 



The parameter uq can be determined from the condi- 
tion ( [4.4| ) which is equivalent to the requirement D{tp) = 
0. The main contribution to the integral comes from the 
region where h'{s) ~ a{s) = which is again the case for 
s ~ ta- With n{s) « ng the result ( 4.24 ) is recovered 
and HQ indeed agrees with the value derived within the 
QFDT-solution. 

It is remarkable that the longest time scale in phase 
A is not the external time scale v^^ but rather t'a ~ 
v^^^^', which is longer. The ultimate reason for that is 
the behavior of K{t) = ki < at the border between 
the plateau and the asy mpto tic regime. This value also 
enters the exponent 77, ( 4.74 ). 



4.6.2 Phase B 

This is no longer the case in phase B. The choice 

ta - y-^ 



(4.83) 



allows to rewrite the complete set of mean field equations 
for t ^ ta ya terms of functions of t/ta only, especially 



n{t) = n{t/ta) . 



(4.84) 



In the plateau region n{t) is given by (4.5£) and the 
dependence on c drops out for cla{t/tx) ~ 7r/2, or with 



(|4.64|) for t - e^'^Hp. Matching with ([4. 84]) yields 



ta^e^'^Hp. (4.85) 



On the other hand, ( 4.69| ) results from matching q{t) 
and therefore 

This means c = c{v) — > for v ^ as proposed earlier. 
With ( 4.64 ) the intermediate time scale t^ is 



t, 



(1-0/2 ,,-(1-1-0/2 



(4.87) 



The force necessary to sustain the velocity v has to 



be calculated from (2.16). The main contribution comes 
as before from the region where 1 + i>{s) — a(s) = 0, 
which is the case for s ^ ta- Using the scaling form 



q{t) = q{t/ta), dp), and (4.84) yields 



F^Fr, 



= -2/3 / da; x f" {q{x) + x"^) 



1 -I- n{x) 



(4.88) 



This does not depend on v and therefore a finite pinning 
force Fp exists, which has to be overcome in order to set 
the particle in motion. 
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4.6.3 Phase C 

In phase C again 



(4- 



has to be chosen and matching to ( 4.66| ) requires = ta- 
Otherwise the same arguments as above (phase B) hold 



and there is again a finite pinning force given by (4.8 
The main difference between phase B and C is the value 
of n{t) for tp <^ t <^ ta, which is finite in phase B and 
essentially zero in phase C. 

It is again of interest to compare this result with the 
one-dimensional case [|l|, ^ where a finite pinning force 
is found for 7=5, whereas the pinning force diverges 
for 7 < i. This divergence is due to the fact that no 
short distance cutoff in the correlation of the disorder 
is used in the one-dimensional calculation. This is also 
the reason why the phase boundaries in this case do not 
depend on temperature. 

4.7 The phase diagram 

Let me summarize the results regarding the phase dia- 
gram and the dependence of v{F). Depending on tem- 
perature T and exponent 7 several phases have been 
found. The above considerations yield the phase dia- 
gram shown in Fig.|] and discussed in more detail be- 
low. Different dependencies of the velocity on the driv- 
ing force are observed in different phases. Examples are 
shown in Fig.n3. 



10^ 



10-3 



10-' 
10-1 



\v r=1.5 (0.75) 






0.40 : 

0A2 
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Figure 15: Drift velocity as a function of the applied 
force in phase O, A, B and at the boundary between 
phase O and A. 

The curve marked .75/. 7 belongs to 7 — 0.75 and 
r = 0.7, which is in phase B. This point is also used for 
the examples presented in Section 3.2. Pinning is clearly 
observed. With increasing force v{F) becomes linear in 
F indicating a finite mobility m — v/F. Its value is 
determined by the first term in ( 2.16| ), which means that 
the influence of the random potential vanishes for high 
velocities. 



The curves marked 0.40 and 0.42 correspond to 7 = 
1.5 and T = 0.4 and T = 0.42, respectively. Both points 
are in phase A. The temperature T — 0.4 is the second 
example in Section 3.2. A power law dependence of v(F) 
i s fou nd in accordance with the analytic investigations 
(|4.8l|) . For 7 = 1.5 and T = 0.4 one finds r/ = 0.053. 
With increasing force again the free mobility is found. 
The curve marked 0.44 is right at the transition line 
7 = 1.5, T = Tc{l.5) = 0.4387, and rj = 0.19 is found. 
The remaining curves marked 0.46 and 0.48 are in the 
high temperature phase O with T = 0.46 and T - 0.48, 
respectively. For v —^ a finite mobility is found, which 
is much smaller than the free value. It vanishes at the 
critical temperature Tc- 

4.7.1 Drift phase O 

For T > Tci'y), ( |4.14| ), and 7 > 1 a finite mobil ity m = 
v/F is found for ?; ^ 0. It is given by ( 4.16 ) and it 
vanishes at the boundaries of this phase. There are no 
long time scales. 

4.7.2 Creep phase A 



This phase is characterized by uq < given in (4.24) 
and ( [1.56 ). Its boundary with phase O is T = Tc(7) 
and the boundary to phase B is determined by ng = Ux- 
This yields with (U), ( p^ ) and ( |456| ) 



f'iqc)f"iqc 



2(l + i^.). 



(4.90) 



Inserting ( p.5[ ) one finds that this phase boundary is de- 
termined by 7 = 7c, with 



1 



7c = 



1.3044. 



(4.91) 



Fo r ?; — > one observes creep in the form v ~ F^^^ , 
(|l8l|) , with 77 for 7 ^ j^- 

4.7.3 Pinning phase B and C 

This phase exists for 7 < 1 or T < Tc and 7 < 7c. It 
has a finite pinning force which has to be overcome in 
order to set the particle in motion. Phase B and C are 
distinct only by their value of n{t) for tp <^ t <^ ta- In 
phase C one finds n{t) Ri for t <^ta which means that 
the FDT holds up to this value, whereas it starts to be 
violated already for t ^ ta phase A and B. 

5 Discussion 

The first aspect of this investigation deals with the mo- 
tion of a particle in a correlated random potential with 
power law decay of the correlations under the infiuence 
of an applied driving force. The phase diagram shows 
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a phase with finite mobihty, a creep phase and pinning 
phases. Similar behavior is found in a one-dimensional 
model [|l|, H indicating that the transitions found are 
not an artifact of the mean field treatment, which be- 
comes exact in the limit of infinite dimensionality. In 
the creep phase v{F) obeys a power law. Such a behav- 
ior is also found in the one-dimensional case, but only 
at the boundary between the drift and the creep phase, 
which is otherwise ruled by a stretched exponential law. 
The pinning phase in the present case has a finite pin- 
ning force, which is also found in the one-dimensional 
case, but only at the boundary between creep and pin- 
ning phase. Otherwise the pinning force diverges in this 
calculation, which can be traced back to the absence of 
a short distance cutoff of the power law decay of the 
correlations. 

The numerical results indicate the existence of sev- 
eral scaling regimes which are verified by analytic in- 
vestigations of the asymptotic properties in the limit of 
small drift velocity. These regimes are, with increas- 
ing time, the FDT-regime describing a local equilibrium 
within one of the valleys of the energy landscape, an in- 
termediate plateau regime where the correlation function 
q{t) stays close to the EA-order parameter qc and where 
the characteristic time scale tp{v) ~ v^, and the asymp- 
totic regime with characteristic time scale ta ^ v^^. In 
the creep phase A this asymptotic regime is ruled by two 
time scales, ta ~ v~^^^ and t'^ ~ v~^~^. The exponent 
rj also determines the power law of v{F). 

The numerical calculations have been performed for 
a wide range of velocities including values as small as 
10~^° and over times ranging from 10^^ to lO'^^. This is 
necessary in order to deduce the full asymptotic behavior 
and even at these extreme values part of the structure 
is not yet fully developed. 

The second aspect relates to glassy non-eqilibrium 
dynamics of mean field models. The model investigated 
here has several advantages in this respect. It uses Lan- 
gevin dynamics which is certainly easier to handle than 
for instance Glauber dynamics in Ising type models. De- 
pending on 7, continuous as well as discontinuous ergod- 
icity breaking transitions arc found. With applied exter- 
nal force a stationary non-equilibrium state is reached 
where correlation and response functions depend only 
on time differences. The inverse velocity plays the 
role of an external long time scale. 

The replica treatment of this model |Q, ^ predicts 
transitions between a phase with continuous replica sym- 
metry breaking for 7 < 7c = 1, a IRSB- phase for 7 > 1 
and T < Tc^msB, and a replica symmetric phase for 
7 > 1 and T > Tc^msB- The present phase diagram 
differs in the sense that Tc > Tc^msB and 7c > 1. A 
difference in Tc has been observed before in models with 
discontinuous transitions ||, |l3|, |l^, |l^ . This was traced 
back to the fact that the states contributing most to the 



static replica calculation are different from those relevant 
for dynamics and are not accessible within finite time in 
the thermodynamic limit. The same now appears to be 
true for continuous transitions as well. 

There have been several proposals regarding the long 
time dynamics. For the SK-model Sompolinsky and Zip- 
pelius ^ proposed a hierarchical or ultrametric or- 
ganization of long time scales. This can be rephrased 
as the postulate ||l^ that correlation and response func- 
tions can be expressed as functions oi x{t) = 1 — Ini/lnt 
where t is some long external time scale. It has, however, 
been shown that this leads to inconsistencies . In the 
present formulation this requirement means 1^(1:) — for 
times where the hierarchy exists. This is not observed. 

The assumption of an ultrametric hierarchy of time 
scales leads to results which are identical to those ob- 
tained by Parisi's continuous replica symmetry breaking 
scheme ||l^. One of the quantities to be compared is 
the probability of overlaps P{q) which is the derivative 
P{q) = of a fimction x{q) which in dynamics is 

given by 

Within this scheme only x{q) is determined, whereas cor- 
relation and response functions are not unique for long 
times. This is in contrast to the present investigation 
where correlation and response functions are unique for 
all times. 

The function x{q) obtained from replica theory and 
from the present investigation are compared in Fig.p^. 
The difference between the results shows again that dif- 
ferent states are of relevance in dynamics and replica 
theory. 




Figure 16: Function x{q) in phase B, see text. The 
results obtained in the present investigation (solid) and 
from replica theory (dashed) are shown. 

For systems with discontinuous transitions ^ |l^, 
[T^ the QFDT-solution has been proposed. It requires 
n{t) uq for t ^ t. This is actually found in phase 
A. Within this scheme nothing can be said about the 
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generation of internal long time scales. Furthermore, the 
resulting correlation and response functions are again 
not uniquely determined, contrary to the present results. 

There has been considerable interest on mean field 
dynamics of spin glasses and related systems with non- 
equihbrium initial conditions ^ |l], |2[j2|, ||, |2|, 
In the work of Cugliandolo and Dean [p5|| an exact so- 
lution of the spherical SK-model is reported. It shows 
aging phenomena although a replica calculation of this 
model does not require replica symmetry breaking. The 
remaining papers deal with models which have phases 
with broken replica symmetry (1-step or continuous). 
Various proposals are made for the long time proper- 
ties of correlation and response functions depending on 
two time arguments t and t' . Again for i — t' — > oo 
and t' ^ oo with finite {t — t')/t' the resulting corre- 
lation and response functions could not be determined 
uniquely. This is likely to be due to an incomplete analy- 
sis of the plateau regime where t—t' ^ t''' with < C < 1- 

A comparison with the present results certainly has 
to be taken with care because the specific non-equilibri- 
um situation is different. Nevertheless using t' as exter- 
nal time scale one expects that at least the properties 
for times t — t'<^t' can be compared. This includes the 
plateau regime. A careful investigation of this regime 
and the kind of FDT- violation taking place there seems 
necessary in this case, too and this is likely to remove 
the arbitrariness in the correlation and response func- 
tions found so far. 



A Dynamic mean field equations 

Time dependent expectation values of products of the 
variables q and q can be represented as path integrals 

{0{q, q)) = Jv{q,q}0{q, q) 

xexp(-^{5°+5;i}) (A.l) 

a 

with 

SZ^ f dt(e„(t)2-|-ie„(i)[(9t-H^o)ea(t) 

(A.2) 



and 



dV{eit)) 



dga{t) 

Averaging over the disorder with (|2.4[) results in 



(A.3) 



Ecu , q _ _ 1 Ov Qv 

^ Dv — 2 '^a'^iB 



(A.4) 



and the averaged action is 

=0' I dt dt' [f'{x{t,t')) Qo.{t)Q^{t') 



(A.5) 



with 



Ql{t) ~ Qiit')] 



(A.6) 



This allows to calculate correlation functions ( |2.6| ), re- 
sponse functions (2.7), using 



and other expectation values like W{t), given by (2.4) 
and (|t^). 

In mean field theory, see e.g. ^ , an effective action is 
introduced replacing appropriate terms in the averaged 
action by their expectation values such that the whole 
problem separates. This effective action is in the present 
case 

S^ff = S° + ifdt dt' W{t - t')Q^{t) go.it') 



(A.; 



+ ij dtj dt' w{t-t')r{t-t')ga{t) 
X [gc.{t) - ga{t') + VNv{t - t')(5„i] . 
The identity 



/.^f .^ SS{g, g) ^ , SO{g, g) ^ 



6m 



(A.9) 



can be used to derive the equations of motion and the 
driving force. The choice 0{g, g) = g{0) and the require- 
ment of causality r{t) =0 for t < yields for a > 1 the 
first mean field equation (2.8). With 0{g,g) = g{0) 
for a > 1 and (|]^) the second mean field equation 
(2.9) is found. The driving force (2.16) is obtained from 
0{g, ^) = 1 for a = 1. 
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